Nonlinear a Model Method for the J\-J% Heisenberg Model: 
Disordered Ground State with Plaquette Symmetry 
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A novel nonlinear a model method is proposed for the two-dimensional J1-J2 model, which is 
extended to include plaquette-type distortion. The nonlinear a model is properly derived without 
spoiling the original spin degrees of freedom. The method shows that a single disordered phase 
continuously extends from a frustrated uniform regime to an unfrustrated distorted regime. By the 
continuity and Oshikawa's commensurability condition, the disordered ground states for the uniform 
J1-J2 model are plaquette states with fourfold degeneracy. 
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The two-dimensional (2D) J1-J2 model is a frustrated 
Heisenberg model with nearest neighbor (NN) and next 
nearest neighbor (NNN) antiferromagnetic exchange in- 
teractions on a square lattice. The model with spin mag- 



nitude S 



is realized in mother materials of cuprate 



superconductors, La 2 CuO_4, YBa 2 Cu0 6 and Sr 2 Cu0 2 Cl 2 
as small- J 2 systems 0, @- Recently found materials, 
Li 2 VOSiC>4 and Li 2 VOGeC>4, are also described by the 
model in the case of J 2 /Ji~l 0,0]. A particular inter- 
est for the Ji-J 2 model is in a gapful disordered state, 
which may be formed by frustration under strong quan- 
tum fluctuations 0. The subject has been theoreti- 
cally investigate^ by various methods 0: e. g. spin 



wave theories [a, |7|, |8( , nonlinear a model (NLSM 
ods 0, j numerical diagonalizations [ill [l 



quantum Monte Carlo (QMC) simulations |1.1ll( 
ries expansions 0,0, H3' an d variational methods [2l|. 

For a system only with the NN interactions (J 2 =0), 
the ground state is believed to have an antiferromag- 
netic (AF) order. The NNN exchange interactions are ex- 
pected to induce strong frustration to break the AF order 
and to form a disordered ground state around J 2 / Ji = 
0.5. A current leading QMC calculation 0, sup- 
ports the disordered phase with spin-gap for J 2 /Ji>0.4. 
Accepting this result, the issue is the character of the 
ground state in the disordered phase. Candidates ex- 
amined in recent several years are a uniform resonating- 
valence-bond (RVB) state |2lj |. a plaquette state 
a dimer state 0, 113, 0] , and a state both with dimer 
and plaquette structures ^t|> their degeneracies are 1, 
4, 4 and 8, respectively. Although Oshikawa's commen- 
surability condition |22j is useful to restrict possibilities, 
it does not completely select one; e. g. it requires that 
a uniform RVB ground state with spin-gap is accompa- 
nied with gapless singlet excitations. The character of 
the ground state is still under debate. 

A disordered state is formed also by distortion in 
the exchange constants, even if there is no frustration 
(J 2 = 0). For a plaquette-type distortion, a disordered 



state interpreted as a 2D array of plaquette-singlets is 
formed |23j. Here it is a question whether the disor- 
dered state by frustration is essentially the same as that 
by plaquette-type distortion. If it is the same, a disor- 
dered phase continuously extends from a regime of strong 
frustration and weak distortion to a regime of weak frus- 
tration and strong distortion in a parameter space. How- 
ever, if not, there exists a phase boundary between them; 
then the ground state of the uniform J\- J 2 model is not 
plaquette-like. Hereafter we consider the Ji-J 2 model 
which is extended to include a plaquette-type distortion. 

Among various methods to analyze spin systems, an 
NLSM method is effective to clarify their characters. The 
first successful example appeared in one dimension. A 
uniform spin chain with NN interactions is mapped onto 
an NLSM with an appropriate topological term |3| ■ In- 
homogeneous spin chains with periodicity are treated by 
refined and extended NLSM methods |2l|2^|. For 2D 
systems, an NLSM without topological term is derived for 
J 2 = 27] . For J 2 7^ 0, Chakravarty et al. analyzed 
2D NLSM which represents the uniform J\- J 2 model. By 
applying a renormalization group (RG) method to the 
NLSM, they constructed a standard theory for the quan- 
tum phase transition. 

Despite the success, there remains ambiguity in the 
correspondence of a derived NLSM to the J\- J 2 model. 
If one use a naive mapping in literature, a single spin 
variable is replaced by the sum of two new variables rep- 
resenting a slowly varying AF motion and a rapid fluctua- 
tion. This is not justified because the number of indepen- 
dent variables is abruptly increased. Although the map- 
ping may phenomenologically produce the correct NLSM, 
there is no way to confirm the correctness within the 
NLSM method itself. Further the increase of the degrees 
of freedom leaves ambiguity for the choice of the cutoff. 
In one dimension, the problem of the degrees of freedom 
has been overcome in generalized formulations [25L l2( : il| . 
However such a reasonable theory in two dimensions has 
not been proposed. To construct a qualified 2D NLSM 
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method is a purpose of this Letter. Using the NLSM 
method, to determine the character of the ground state 
for the 2D J\-J% model is the final purpose. 

The J1-J2 model with plaquette-type distortion is rep- 
resented by the Hamiltonian: 



(a) 



H 



E 



J: 



2:ik 



S,; ■ Si 



(1) 



where is the spin of magnitude S at site i. The first 
and the second summations are taken over NN and NNN 
pairs, respectively, in a square lattice. Ji-ij takes J\ or 
J{, and J 2; ifc does J 2l J 2 or J 2 as shown in Fig. [Ja). 
The system is reduced to the uniform J1-J2 model when 
Ji = J[ and J 2 = J 2 = J 2 ■ In the limit of J[ — J 2 — 
J 2 = 0, the lattice is an assembly of isolated plaquettes 
each of which consists of four spins connected by J\ and 
J2 (Fig. Gib)). Also, in the limit of J x = J 2 = J' 2 ' = 
0, the lattice is an assembly of another kind of isolated 
plaquettes; each consists of four spins connected by J[ 
and J 2 (Fig. Efc)). Hamiltonian JJJ is invariant under 
the simultaneous exchanges of Ji and J{, and of J 2 and 
J 2 . The symmetric case of J\ = J{ and J 2 — J 2 includes 
the uniform J\-J 2 model. 

We consider the quantum Hamiltonian in the clas- 
sical Neel ordered region. The expectation value of Sj 
for a spin coherent state at imaginary time r is given as 



-l) j Siij(T) with n 2 = 1, 



(2) 



where (— 1) J is a symbol taking + or — depending on the 
sublattice which the jth site belongs to. The partition 
function is then written in a path integral formula as 

Z=(D[n j {r)]Y[5{^{T)-l)e- A . (3) 
J j 

The action A at temperature 1/(3 is given by 

rP 

A = t£y7-l) J 'u7[n J ]+ / (1tH(t). (4) 
Jo 



The first term is the Berry phase term with the solid 
angle u>[nj] which the unit vector iXj(t) forms in period 
(3. H(t) in the second term is given by 



H <r) = ^E^Mr)-^)] 2 

- \s 2 J2 J2; ik [Mr)-n k (r)} 2 , (5) 

«i,fc» 

where the constraint n 2 (r) = 1 in the <5-function of 
Eq. J2J has been used. Hereafter we do not explicitly 
write the r dependence of (r) . 

We adopt a plaquette of Fig. Qlb) as a unit of transfor- 
mation, and call it a block: we would choose another kind 
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FIG. 1: (a) Lattice of the J1-J2 model. Lattice sites are 
denoted by small circles and exchange constants are by various 
lines between circles, (b) A plaquette consisting of four sites 
connected by Ji and J2; this is a block which is a unit in the 
NLSM formulation. A variable nj in the pth block is relabeled 
as n ,u/ (p), where n and v take + or — . The value of (/x, v) at 
each site is shown, (c) Another kind of plaquette consisting 
of four sites connected by J[ and J 2 ■ 



of plaquette in Fig.^c) as a block. We relabel four vari- 
ables, rij's, in the pth block as n ++ (p), n H (p), n~ + (p) 
and n (p) , as shown in Fig. ^b) . By analogy with the 
one-dimensional case |26j |. we transform them as 

n"» = m(p) + a[fivL (p) + M Li(p) + i/L 2 (p)]. (6) 

Here Lo(p), Li(p) and L 2 (p) describe small fluctuations 
around m(p). According to the variable transforma- 
tion, four original constraints, [n MI/ (p)] 2 = 1 (/i, v = 
±), are changed to four new constraints, m 2 (p)=l and 
m(p)-Lg(p) = (q = 0, 1, 2). Thus we obtained a new set 
of variables, the number of which is the same as that of 
the original variables. This plaquette-based transforma- 
tion is inevitable to keep the original degrees of freedom 
even in the uniform J\-J 2 model. 

In the continuum limit, the first term of the 
action J3J is written as iS J2 p v p,vw\a. llu {p)\ = 
i(S/a) JdrcPr Lo{mx<9 r m) with lattice spacing a. For the 
second term of Eq. (gj, we substitute Eq. © into Eq. (J5J 
and take the continuum limit. Thus, to the leading order 
of derivatives and fluctuations, we have the field-theoretic 
action 



A = S 2 J drd 2 r{^-L -(mxa r m) 



+ J [{d x m) 2 + {d y m) 2 - 2d x m ■ U - 2d y m ■ L 2 ] 



2(J 1 + J[)L 2 + (J Q + 4)(L 2 +L 2 )} 



(J) 



with J = Ji - J 2 - J'J, and J' Q = J[- J' 2 - J' 2 \ This 
action includes all the low-energy excitations surviving 
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the continuum approximation, since the original degrees 
of freedom are not spoiled in the variable transformation 
©. In Eq. 0, L , h 1 and L 2 are massive fields |28| . 
so that they are irrelevant to a symmetry change of the 
ground state. 

Now we integrate out the partition function for the ac- 
tion J7J) with respect to massive fields Lo, Li and L 2 . The 
resultant partition function contains the NLSM action: 



A e g = / drd 2 r 



1 



5'- 



1 1 

Jo Jn 



-1 

[(d x m) 2 + (d y m) 2 : 



There appears no topological term even if the NNN in- 
teractions exist. The bare spin wave velocity is read as 
V = 2V2S*a(Ji + J[)^ 2 (1/J Q + 1/J^y 1 / 2 . Action A cS 
keeps the original invariance against the simultaneous ex- 
changes of Ji and J[ , and of J 2 and J 2 , meaning that the 
same action is obtained if we use a plaquette in Fig.^c), 
instead of Fig.^b), as a block. This result reflects that 
the variable transformation 10 does not restrict the spin 
motion to form a singlet on the plaquette of Fig.QIb). 

We apply the RG analysis by Chakravarty et al. |9j 
to the present NLSM. We first introduce rescaled dimen- 
sionless coordinates, xq = Avt, x\ = Ax and x 2 = Ay, 
with a momentum cutoff A of order a -1 . The NLSM 
action JHJ is then rewritten as 



1 

A c s = 

2 5o 



d 3 j 



dm 

dx u 



(9) 



with coupling constant g = V2AaS~ 1 (Ji + J() 1 / 2 (l/ J + 
1/Jg) 1 / 2 . By RG equations up to one-loop approxima- 
tion, the quantum phase transition from the AF ordered 
(Neel) state to a disordered state takes place at go = Air. 
Rewriting this, the phase boundary in the space of the 
exchange parameters is given by 



(J, (J- • 
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2 
A 



with A 



/ Aa 



\2ttS 



(10) 



Parameter A represents the strength of quantum effect; 
A = in the classical spin limit. 

To make the NLSM method complete, we determine 
the cutoff A by considering the number of degrees of 
freedom for the square lattice. The variable m is orig- 
inally defined for each block of size 2a x 2a (Fig. ^b) 
and Eq. ©), and is taken a continuum limit. Hence the 
correspondence of the momentum spaces is expressed as 



(■7r/a) = 7rA , or the cutoff is given by A = 



Thus 



Eq. (|10|1 unambiguously determines the phase boundary 
between the ordered and the disordered phases. 



In the uniform limit (Ji = J[, J 2 = J 2 = 



the 



system depends only on frustration parameter a = J 2 / Ji 
and Eq. I|ll)fl is reduced to a — i — A. Hence, for S = | 
with A = 



a c ~ 0.18. Thus the NLSM method succeeds in pro- 
ducing a critical value satisfying < a c < ^ without 
any additional assumption or interpretation. The value 
is smaller than ^0.4 estimated by the QMC simulation 
(ljlllij ]. The deviation reflects the difference between the 
dispersions for spin-wave excitations in the lattice and 
the continuum models, and may be reduced by adjusting 
the cutoff. Since we aim at inspecting the continuity of a 
phase, we do not need such a phenomenological adjust- 
ment. 

In the limit of no frustration (J 2 — J' 2 — J 2 = 0), the 
plaquette distortion may cause an order-disorder transi- 
tion. We denote the strength of the distortion by distor- 
tion parameter^ defined as J[ = (1 — 7) J\. Then Eq. i|10|) 
produces the critical value 7 C = 2 — A -1 + \/\~ 2 — 2A _1 . 
This value decreases from 1 to as A increases from to 

_L 

2 ■ 

We now examine the continuity of the ground state 
between both the limits above. To be concrete, we pa- 
rameterize the exchange constants as J[ — (1 — 7) J%, J' 2 
= (1 - 7) 2 J 2 and J' 2 ' = (1 - 7) J 2 for < 7 < 1. Equa- 
tion ijTTJI) for the phase boundary is reduced to a simple 
form as a = (2 — 7) _1 — 5-^(2 — 7)(1 — 7) _1 - The phase 
diagram in the 7-a parameter space is shown in Fig. [21 
The bold line with S—oo is the classical phase bound- 
ary between the Neel and the colinear phases [28|, The 
phase boundary of S=^ between the gapful and the gap- 
less phases for variable m is the thin solid line; the state 
above is gapful, while that below is gapless correspond- 
ing to the Neel (AF) ordered state. Boundaries for other 




0.2 



0.0 



plaquette 



AF order 



S = 1/2 



0.0 



0.2 




1.0 



7r/a, the critical value for a is given as 



FIG. 2: Phase diagram in the space of distortion parameter 
7 and frustration parameter a (=J2/Ji). The bold solid line 
for S=oo separates the classical Neel and the classical colin- 
ear phases. The region between the bold solid and the thin 
solid lines is the gapful plaquette phase for S=^. The phase 
boundaries for S — 1, | and 2 are also shown by dashed lines. 
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spin magnitudes S are also shown by dashed lines. 

The gapful region of m in Fig. [5] extends continuously 
from the uniform limit on the a-axis (7=0) to the limit 
of no frustration on the 7-axis (a=0). Remembering that 
fields Lo, Li and L2 are gapful, there is no gapless excita- 
tion throughout the region whether it is triplet or singlet. 
Hence, the whole gapful region in Fig. [21 is a single dis- 
ordered phase. In particular, the phase continues to the 
point of (7, a) = (1, 0) 29j. Hence a disordered ground 
state on the a-axis finally continues to the ground state 
of the assembly of isolated plaquettes. 

Thus there remain two possibilities for a disordered 
ground state of the uniform J\-Ji model, which is on the 
a-axis in the phase diagram (Fig.EJ). First, the transla- 
tional symmetry may be spontaneously broken; then the 
ground states are fourfold degenerate and one of them 
continues to the ground state at (7, a) — (1, 0). Second, 
the symmetry may not be spontaneously broken; then 
the ground state is unique and is a uniform RVB state 
with strong fluctuations of plaquette-singlets. However, 
the second possibility is excluded by Oshikawa's com- 
mensurability condition |22l |. Applying it to the present 
case, a uniform ground state with triplet excitation gap 
must be accompanied with other gapless excitations like 
singlet ones. Such gapless excitations do not exist as we 
have already shown. We therefore conclude that the dis- 
ordered ground states for the uniform J\-J-i model are 
fourfold degenerate plaquette states with spontaneously 
broken translational invariance. 

Finally, we discuss possible experiments to detect the 
plaquette state for materials with J2/J1 ~ 0.45 which 
will be hopefully found in future. In a realistic layered 
material, the uniform J\-Ji model is accompanied with 
at least weak three dimensionality. Hence, at a finite tem- 
perature, the system will spontaneously break the trans- 
lational symmetry to fall into a plaquette phase. The 
appearance of a spin-gap at the temperature will be ob- 
served. The characteristics of the plaquette state appear 
in the dispersion relation, which reflects the invariance 
for the translations of 2a in the x- and the y-directions. 
They will be observed in neutron scattering experiments. 
In a material where the spin system weakly interacts with 
the lattice, the spontaneous symmetry breaking for the 
spin degrees of freedom may induce a plaquette-type lat- 
tice distortion, which will be observed by X-ray diffrac- 
tion. Such a distorted system may correspond to a point 
deviated from the a-axis in the plaquette phase of Fig.0 
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